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NUMERICAL APPROXIMATION OF POSITIVE POWER 
CURVATURE FLOW VIA DETERMINISTIC GAMES. 


HEIKO KRONER* 

Abstract. We approximate the level set solution for the motion of an embedded closed curve 
in the plane with normal speed max(0, where k , is the curvature of the curve and 2 < -y < i by 
the value functions of a family of deterministic two person games. We show convergence of the value 
functions to the viscosity solution of the level set equation and propose a numerical scheme for the 
calculation of the value function. We illustrate the convergence properties of the scheme for different 
parameter values in some example cases. 


1. Introduction. In [24] R. Kohn and S. Serfaty present a game theoretic 
approach to the curve-shortening flow, i.e. the mean curvature flow of curves in 
the plane. The (time-dependent) level set formulation of this flow is a degenerate 
parabolic, possibly singular equation and is interpreted as the limit of the value func¬ 
tions of a family of discrete-time, two-person games. It is proved that these value 
functions converge to the viscosity solution of the level set equation, cf. [24l Theo¬ 
rem 2]. Furthermore, using a ’minimal exit time game’ a rate for the convergence 
of the value functions of this game to the time-independent level set formulation of 
the (mean) curvature flow is proved in |24] . We note, that the corresponding conver¬ 
gence result (without rates) for the positive curvature flow (here the normal speed 
is the curvature if it is nonnegative and zero otherwise) is given in [241 Theorem 5]. 
In contrast to the curve shortening flow not much is known about the behavior of 
a non convex initial curve in the plane moving by positive curvature flow. But the 
latter has been used for image processing, see the remarks in [531 page 16] and [55] . 
In form of examples modified flows are presented in [241 Section 1.7] which have a 
game theoretic interpretation. But proofs of convergence for the value functions of 
the corresponding games for these modified examples are not given therein. Among 
these is a game, cf. [54| Example 3, Section 1.7], which can handle the case of the 
curvature flow of a convex curve in the plane at which the normal speed is given - 
instead of the curvature - by a power 7 > 0 of the (mean) curvature. We call this 
power curvature flow (PCF) (and the analogous flow for hypersurfaces power mean 
curvature flow (PMCF)). PCF is only a well-defined flow when the initial curve is 
convex. In [53] five open questions are formulated, cf. [241 Section 1.8]; among these 
is the question of numerical advantages or disadvantages of this approach. 

In this paper we consider for ^ < 7 < 1 the more general, so-called positive 
power curvature flow PPCF for which the normal speed is given by (where k+ = 
max(0, k)) and which coincides with PCF if the initial curve is convex. Our aim 
is to approximate this flow by a family of value functions of suitable games. We 
show convergence of the value functions to the unique viscosity solution of the time- 
dependent level set formulation of PPCF and present a numerical algorithm for this 
approach. We demonstrate the convergence of the algorithm for some example cases 
and analyze the effect of different choices for the numerical parameters. 

Let us relate the content of our paper to [ 53 ]. There is no game in [ 53 ] for PPCF. 
The game in [531 Example 3, Section 1.7] is suggested therein to approximate PCF 
in case 0 < 7 < 1. The corresponding limit equation [531 (1-23)] has the (formal) 
disadvantage that its left-hand side becomes — c» if there are non convex zero-level 
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sets (in case they are curves) of the continuum solution (e.g. at the initial time outside 
of the initial curve). By considering PPCF we avoid this convexity issue, especially 
we can use non convex initial curves. The game we use for the approximation of 
PPCF is essentially the one in Example 3, Section 1.7] but it turns out that it 
is necessary (at least for our proof) to use a slightly different definition of the value 
function. The difference is that we take the minimum in the definition of the value 
function over a (e-dependent) compact subinterval of ( 0 , oo) instead over ( 0 , oo), cf. 
( 1 ^ . This difference is not due to the fact that we consider PPCF instead of PCF but 
is of general nature. The initial and endpoints of the subintervals have to satisfy 
some technical relations, cf. which imply that we need the restriction 

i < 7 < 1. The interesting point is that i is exactly the critical value at which the 
convergence behavior of PCF changes, see Theorem 11.21 and the following remarks. 

Concerning a numerical application of the Kohn-Serfaty approach to the curve 
shortening flow (i.e. 7 = 1 ) we refer to |llj . In addition therein a semi-Lagrangian 
approximation of the curve shortening flow is considered. A semi-Lagrangian approach 
for the affine curvature flow case which corresponds to 7 = | is content of m- 

We give some references for the numerical approximation of flows by powers of 
the curvature (which do not use the Kohn-Serfaty approach). A numerical approx¬ 
imation of the stationary level set solution of PCF in case 7 > 1 is content of our 
two previous papers [ 2 ^ and m- There we approximate the regularized stationary 
level set equation (used in [36] to prove existence of the stationary level set PMCF) 
by finite elements and prove rates for the approximation errors. The stationary level 
set equation of mean curvature flow is also content of nni. 

There are only few numerical papers which deal with PCF for 7 7 ^ 1. 

In [7] and [5] stability bounds for a finite element approximation of PCF for 
curves is derived, see also |9]. In [30] the parametric formulation of the evolution of 
plane curves driven by a nonlinear function of curvature and anisotropy is considered. 
Concerning PCF from the perspective of image processing we refer to mm.m.m 
and 

For estimates of the regularization error for the instationary level set formulation 
of mean curvature flow we refer to |31j and m- For more general references to the 
numerical treatment of motion by curvature we refer to m, [m, [H], [32], [E] , [20], 
[ 33 ] and [40] where the instationary level set formulation is considered. 

Our paper is organized as follows. In the remaining part of the present section 
we formulate some known results about the behavior of curves evolving under PCF. 
In Section]^ we define the value functions u'^. In Section |3| we show convergence of 
the value functions to the solution of the continuous problem. In Section |4| we present 
the algorithm in order to calculate u^. Its application in some example cases as well 
as the effect of the choice of different numerical parameters is content of Section |5| 

As for the positive curvature flow little is known about the behavior of PPCF 
for not convex initial curves. We state some well-known facts about the behavior of 
curves evolving under PCF in case 0 < 7 < 1, see |4] and |3|, and remark that in 
these references also the case 7 > 1 is treated. 

Theorem l.I. Let 7 = 1 and Fq a smooth, embedded closed curve given by an 
embedding xo : —>■ . Then there exists a unigue smooth solution x : x [0, T) —>■ 
M2 of 

(l.I) —X = —K^v (= —Kv) 

where v denotes the outer unit normal, with initial data xq. Ft = x{S^,t) converges 
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to a point p G as t ^ T. The rescaled curves (Ft — p)l •\/2(T — t) converge to the 
unit circle about the origin as t ^ T. 

The assumption of smoothness of the initial curve in Theorem 1 1.1 1 can be relaxed 
to allow Tg to be the boundary of a bounded open convex region and the curves F^ 
are converging to Fg in Flausdorff distance as t —>■ 0, see [5]. 

If 7 = ^ and Fg a smooth, convex, embedded closed curve then Theorem 11.11 
still holds with a different rescaling, namely (Ft — p)/(4(r — t)/3)i converges to an 
ellipse of enclosed area tt centered at the origin. The regularity assumption on the 
initial curve can be weakened to allow boundaries of open bounded convex regions. 
For 7 G [i, 1] we get the following result. 

Theorem 1.2. //7 G 1] o,nd Fg the boundary curve of an open bounded convex 
set in then E2P has a smooth and strictly convex solution for 0 < t < T which 
converges to a point p G R^ and the rescaled curves (Tt — p)/(T — converge 

smoothly to a limit curve F-^ which satisfies k'*' = A (x, v) for some A > 0. 

Curves Fa are called homothetic solutions and in case 7 > ^ all homothetic 
solutions are circles. For 7=5 the only homothetic solutions are ellipses. For 
0 < 7 < i a classification of homothetic solutions is given in [H Theorem 1.5]. 

There are results for non-convex initial curves in and [HJ Chapter 8 ]. 

For 0 < 7 < ^ the isoperimetrical ratio becomes unbounded near the final time 
for generic symmetric initial data, cf. the remarks in [H Section Ij. 

There are results on PMCF in case of n-dimensional hypersurfaces in Euclidean 
space, n > 2, see for example [3S] and m for the parametric formulation and [36] 
for the time-independent level set formulation in case 7 > 1 and |31j (and references 
therein) for the time-dependent level set formulation of PCF in case 0 < 7 < 1. 

Concerning a time-dependent level set equation for PCF we only know the equa¬ 
tion 



in case 7 = 2 n-i > R G N, which is a little bit different from PCF since here negative 
curvatures are allowed. For PPCF the time-dependent level set equation is given by 



for general 0 < 7 < 1. We consider these equations in R^ x (0, 00 ) with initial condition 
m(-, 0) = Mg. We assume that Mg G (^^^(R^)^ |up| < 1 and Mg(a;) = 1 for a; G R^\i3_R(0) 
where i? > 0 is a sufficiently large radius. 

Equations (ESI), (O are fundamental equations in image processing, especially 
in the case 7 = | where these equations are the so-called affine curvature equations, 
cf. [T] and [22) . Equations (HU), (O are included by the general geometric equations 
which are dealt in [H, m and [ 22 ] and therefore for each equation exists a unique 
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viscosity solution. We mention m where convergence with rates of a regularization 
of p.2D is proved. 

2. Game theoretic interpretation and definition of Let be a smooth, 
bounded domain in the plane and T > 0. Let uo be a Lipschitz continuous function 
in the plane with mq < 0 in O, itg > 0 in ]R^\f2 and |uo| < 1. We assume ^ < 7 < 1 
and consider the following continuum equation 

Ut — I Vit|(/5(curv(u)) = 0 
u(-,T) =uo 


curv(u) = — div(Vu/| Vu|) 


( 2 . 1 ) 

where 

( 2 . 2 ) 


is the curvature of the level set and 


(2.3) 


ip{k) 


— IkP, k < 0 

0, K > 0. 


We note, that in this sense the curvature of the unit circle is —1 because in {x G 
: |a;| — 1 = 0} we have curv(|a;|) = — div = —n < 0 in view of (12.21) . 

Furthermore, m is not PPCF but after transforming the time variable t ^ T — twe 
get the correct equation, i.e. PPCF (11.31) . As already explained above our equation 
(EH) differs from the continuum equation in [241 Section 1.7, Example 3] in the fact 
that the latter requires Ut = 00 if curv(M) > 0. 

The goal is to approximate m by the value functions of certain deterministic, 
time-discrete two person games. For the convergence proof it will be crucial that for 
K < 0 we can write 

(2.4) (p{k) = sup ( ]-KS^ - f{s)) 

s>0 / 

where 

(2.5) f{s) = c^s^, = (1-7)(27)t^. 


The value functions e > 0, are defined for x G fl and 0 < t < T by 

(2.6) u’^{x,t) = min max u'^ix + besv + (3e'^ f(s)v'^ ,t + e'^) 

lhll=i.se/. b/3e{-i.i} 


where 


(2.7) 
with 

( 2 . 8 ) 
and 
(2.9) 


— [^ 5 ' 


1-7 ■ f ^ 

—— < Oi < mm 1 1 

27 


, 27 ^ 1 

0 < a2 < min ai- -1, - 

1 — 7 3 
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The difference from our definition to the definition in [531 Example 3, Section 1.7] is 
that we minimize over s G instead of s > 0. The first inequality in (12.81) implies 
that the RHS of (12.91) is positive so that there exists a corresponding a 2 ■ The ai < 1 
part of the second inequality in (12.81) is needed for Case (1) in part (ii) of the proof 
of Lemma 13.31 The other part of the second inequality in (12.81) implies that the step 
size is at most a fixed positive power of e. The first part of the second inequality in 
dlSl) is needed for part (i) of the proof of Lemma 13.31 while the second part is used 
for Case (2ba) of part (ii) of the proof of Lemma 15751 Inequality (12.81) implies T' > i 
which is the reason why we have only a convergence result for this case. 

We can calculate numerically by setting = uq and ’playing the game 

backwards in time’, i.e. we calculate u’^(-,T — e^), etc.. 

3. Convergence of the value fnnctions. In this section we show that the 
value functions converge to a solution of m- For completeness we give the following 
definition. 

Definition 3.1. (i) A lower semicontinuous function u is a viscosity supersolu¬ 
tion of if whenever (j){x,t) is smooth and u — (p has a local minimum at (a;o,to) 
we have 

(3.1) eft - |V(()|(/9(curv((())) < 0 
at {xo,to) ifS/(j){xo,to) 0, and 

(3.2) 

at (xo^to) ifV<p{xo,to) = 0. 

(li) An upper semicontinuous function u is a viscosity subsolution of S2.1\) if 
whenever (j){x,t) is smooth and u — 4> has a local maximum at (ccojto) 'xe have 

(3.3) (ft - |V(()|(p(curv(0)) > 0 
at {xo,to) ifV(j){xo,to) ^0, and 

(3.4) Pt>0 

at (xo,to) ifV(j){xo,to) =0. 

We define 


(3.5) 


u{X) = limsup u^iY) 
y->x,e->-o 

uiX) = liminf u'^(Y). 

Y^X,e^0 


u is upper semicontinuous, u is lower semicontinuous and u < u We will show that 
It is a viscosity subsolution of m and that u is a viscosity supersolution of (EH). 
Then from a comparison principle we get that u = u is a solution and that u'^ —^ u. 
The proof strategy starts with an argumentation adapted from |24l Section 4.3] and 
then we distinguish several cases to handle our flow. 

Lemivia 3.2. u is a viscosity subsolution of 

Proof We argue by contradiction. If not, then there is a smooth p such that 
(xo, to) is a local maximum oi u — (p with 

(3.6) (pt - |V(()|(/9(curv((/))) < 00 < 0 at {xo,to) 
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if V(f){xo, to) 7 ^ 0 or we have S/(j){xo,to) = 0 and 
(3.7) ^t<Oo<0 


in (xo^to)- Adding to cj) a nonnegative function whose derivatives at (a;o,to) are all 
zero up to second order, we can assume that (xq, to) is a strict local maximum of u —^ 
in a ^-neighborhood of {xq, to). Let (a;°, t°) —>■ {xo, to) such that u'^(x°, t°) ^ u{xo,to). 

(i) Let us consider the case S/(j){xo,to) ^ 0. We may assume that V(^ 7 ^ 0 in the 
above (5-neighborhood. We construct the following sequence 


(3.8) 




xt = 


=xi 


boeso 


hkf-Sk 


V^0(^o) 

\ mx ^)\ 

|V<^(A^)| 


Poe'^fiso) 


■ f{sk) 


\ mxs )\ 

V4>{xi) 


t^,+e‘ 


|V(/>(X|)|’ 


where the bk and fik are maximizing u'^{X^_^_^) for given Sk and Sk G le is chosen as 
will become clear in the following. Then u^{Xf.) < tt'^(A^_|_j^) in view of (12.61) and 
hence 


(3.9) u^X^o) < ^^{Xl). 

We let X{s) be the continuous path that affinely interpolates between these points 

i.e. X{t) = Xl + ){XIj^i — XI) for -|- ke^ <t<t'l + {k + l)e^, and write 

{x(t),t) = X{t). We have 


(3.10) 

and since 

(3.11) 


m 


4>{X{t)) =dtcl){x{t),t) + X(j){x(t), t) ■ dtx{t) 


_a 

Wt 


Vcj){x{t),t) =X(j)t{x{t),t) + D'^ip{x{t),t)dtx{t) 


dtx(t) = 


bkSk V^0(W|) 

e \mxi)\ 


■/dkfisk) 


mxi ) 

\ mxi )\ 


is constant we get 

(3.12) ^(b{Xit)) = dU{X{t)) + 2V^t{X{t)) ■ dtx{t) + {D^<j){X{t))dtxit), dtxit)). 
Taylor expansion of t (f){X{t)) at -I- ke^ gives 
<l>{Xl^,)-cj,{Xl) 

=£2 {dtdpixD + v^{xi)dtx{t)) 


(3.13) 


+ - {dU{Xl) + 2VMXl)dtx{t) + {DmXl)dtx{t), dtx{t))) 

+ o{e^) 

=e^ (at<).(A^) + |V(/.(X^)|/3fc/(sfc)) 

=6^ (^dt(j){Xl) + \V(j){Xl)\{Pkf{sk) - ^ curv(()5)|xe)^ -f o{€^) 
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in view of 


(3^14) - cutv(^)|D0| = IV0I cliv = (cV^. ^ 

We use the notation 

(3.15) (^s(k) = Ky -/(s) 
for K G R and s > 0 and we get 

(3.16) </)(X|+i) - ldt</,(X^) - |V</)(X|)|(^,,(curv(</))|;,.)} + o(e^) 

We want to show by using (EH) that 

(3.17) dt</>(X^) - |V^(X^)|<p.,(curv(0)|x^) < | 

for small e. Let us assume for a moment that this is already shown then we have 

(3-18) </.(X^+,)-</'(^D<e^y 

and hence 

(3.19) </>(X^) - <P(X^) < < 0. 

Adding this inequality to (13.91) gives 

(3.20) u^(XS) - c/>(X^) < u^iXl) - <j^{Xl) + 

For each e we consider the finite sequence (A^)o</c</c(e) where fc(e) is as follows. Let 

(3.21) U = Bs{xo,to)\Bsi{xo,to), 0 < Si < 6, 
we may assume that Xq G Bs^ (xq, to) and that 

(3.22) m"(Ao) - </>(X^) > u(xo,to) - S)(xo,to) - 
Let 

(3.23) k{e) = minj/c G N : X^ G U}. 

Note, that fc(e) is well-defined since (I3.20|) and (13.221) imply that the sequence (X^) 
leaves Bsi{xo,to) at the latest for k = [e~^]. Hence (for small e compared to S) there 
must be an element X^ &U. 

For a subsequence —>■ {x',t') ^ (a;o,to) and by (13.201) we get 

(3.24) {u-4>){xQ,tQ) < (u-4)){x',t') 

in contradiction to the fact that (a;o,to) is the unique maximum in Bs{xo,tQ). 
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We show (13.171) by distinguishing cases, thereby we assume S to be sufficiently 
small. We set Kq = curv(^)|Xo and = curv((/))|x'. 

Case (1): If kq < 0 then there exists cq = eo(Ko) > 0 so that we can choose 
Sfc € le with 

(3.25) <Psk{Kk) = <f{Kk) 

for all 0 < e < ep. Then we use (ICT . 

Case (2): We assume kq > 0. For small e we have 

(3.26) (ko) < 0 < ip^-o.2 (ko) 

so that there is a Sfc G C (actually not depending on k) with 

(3.27) (^,Jko)=0. 

It follows that (for <5 —>• 0) 

(3.28) tfsk {i^k) -)> 0 = ip{Ko). 

Case (3): We assume kq = 0. 

Case (3a): Let Kk > 0. For small e we have 

(3.29) (Kfc) < 0. 

If 

(3.30) (p^-c2{Kk) > 0 
we choose Sk G le with 


(3.31) (ps^{Kk)=0 

and are ready, otherwise we set Sk = e““^ and get (as e ^ 0) 

(3.32) 0^-/(e-“^) <(^«,(/vfe) <0 


which also implies the claim. 
Case (3b): Let Kk < 0. Let 


(3.33) 


2a2 

1-7 


> fl> 202 


and distinguish the following cases. 

Case (3ba): If \Kk\ < we set Sk = and can estimate 

(3.34) - ^,,{Kk) < + ce“=^. 

Case (3bb): Let \Kk\ > e^. We calculate the maximizing s > 0 of (ps{Kk). For 
general k < 0 the maximizing s > 0 of (/?s (^) is given by 


(3.35) 


= 0 
as 





hence 


(3.36) 



Evaluated for k = Kk we get the maximizing 


(3.37) 


Sk 


f 7-1 

I 27C7 


Kk 


<cei^- 


-y-1 

2 


The so defined Sk lies in if 

(3.38) ^ 
which is the case in view of (j3.33p . 

(ii) Let us consider the case V0(a;o,to) = 0. We define the XI as before by (13.8|) 
but we replace in equation (13.81) the directions and by and Vk 

respectively, Vk an unit vector with 

(3.39) («YV</>(X^))=0. 

The Sk are chosen as will become clear in the following and the /3fc and bk are chosen as 
before, i.e. so that they maximize it'^(X^_|_]^) for given Sk,Vk- Then we have u'^{Xq) < 
u'^{X^) as before. Looking at (I3.13P and (13.111) and adapting it to the present situation 
gives 


cj,{Xl^,)-cl,{Xl) 

( 3 . 40 ) (dtHXD + {yct>{Xl),Vk) Pkfisk)) 

+ ^blsl {D^(j}{Xl)vi,vi) + o{e^). 

Our goal is to ensure that (13.181) also holds in this case. We set 

(3.41) Xl:=\{Vcl,{Xl),Vk)\^0 

as d —>■ 0. If = 0 then we are ready by setting Sk = ■ Let us assume that A| ^ 0. 

We get the following inequality 

(3.42) </)(X^+i) - <).(X^) < {dtcl^iXl) - (Al)} + o(e^) 

where 

(3.43) Al = -Xr^\{D^cl,{Xl)vi,vi) |. 

Case (1): For those which are contained in some sufficiently large interval 
[a, b] C (—oo, 0) (where a, b do not depend on e and k) we choose Sk € h so that 
‘PsdAl) = ip{A%) for small e. 

Case (2): If > 0 we set Sk = e““^ and are ready. 

Case (3): If 0 > > 6 we choose /r according to (13.331) and distinguish cases. 
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Case (3a): If then we set s^ = e and have 

(3.44) 

Case (3b): If |A|A|| > e'^ then |A|| > We calculate the maximizing s > 0 of 
(/3 s(A|) according to p.36l) and get 


(3.45) 


Sfe 


7-1 

27C7 



7-1 

2 


<ce^^A!, 


The so defined Sk lies in if 

(3.46) ^ '^2^ ^ 

which is the case in view of p.33l) . 

Case (4): If < a we choose 


(3.47) 


2Qfl 

1-7 


> M > 7 


2ai 

1-7 


and distinguish cases. 

Case (4a): We assume X% < e^. We set = e“T calculate 


(3.48) 


Afc/(sfe) < ce 


/^+“i 


0 


so that (13.181) follows as well. 

Case (4b): We assume A^ > e^. Then 

(3.49) lA^I < ce-^ 


We calculate the maximizing s > 0 of and get 


(3.50) 


Sk 


f 7-1 
I 27C7 


At 


>ce^- 


7-1 

2 


To ensure that the so defined Sk lies in we need 

(3.51) ai > M— 
which follows from (13.471) . There holds 

(3.52) = ^{AD =-\AIP 
and hence 

(3.53) AI.^.,(7l^) = -Al'-'^l {D^MX^k)vtvt) T ^ 0. 
In all cases we get 

<^(X^.+i) - </.(X^) < 


(3.54) 
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Now, we argument as in (i). □ 

Lemma 3.3. u is a viscosity supersolution of 

Proof. We argue by contradiction. If not, then there is a smooth (j) such that 
(xq, to) is a local minimum of u — cf) with 

(3.55) (ft - |V(/)|v3(curv(0)) > 0o > 0 at (a;o,to) 
if V(f){xo, to) 7 ^ 0 or we have S/(j){xo,to) = 0 and 

(3.56) (ft>Oo>0 


in (xoito)- Again without loss of generality, we can assume the minimum is strict. 
Changing 9o if necessary, we can also find a 5-neighbor hood of (xo, to) in which these 
assertions hold. Again we can find (x°, t°) — {xq, to) such that u'^{x^, t°) u{xo, to). 
Taking sq and the unit-norm vo that achieve the minimum in the characterization 
(12.61) we find 

(3.57) M'^(x°,t°)= max u'^(xo + besoVo +/3e^f(so)vo,to + e^). 

6,/3e{-i,i} 

We set Xq = {x^, t°), Xl = {x + boesoVQ + f[so)vo, to + e^) and inductively 

(3.58) Afc+i = Xl + [bkCSkvi + f{sk)vk,e^) 

where Sk,Vk are chosen recursively (and also depend on e) so that the minimum is 
attained and bk,(3k will be chosen later. We have for all bk,fik S {1, —1} 

(3.59) u^iXl) > u%Xl^,) 
and thus 

(3.60) u^{Xo) > u%Xk). 

On the other hand, extending the X^ into an affine path by affine interpolation as 
before and doing a Taylor expansion gives - adapting (I3.13|) and (13.111) to the present 
situation - 

(/.(A^+i) - (fiXl) = ebkSk {mXl),vi) 

(3.61) + (^dtcfiXl) + (3kf{sk) (V<).(A^), Vk) + ^ {D^cf{Xl)vi,vi)^ 

+ 0{e^). 

Our goal is to show 

(3.62) tf{Xl^,)-f,{Xl)>e^^-^ 
for all sufficiently small (5 > 0. 

(i) Let us consider the case where \/(j){xQ,to) = 0. Let us denote the RHS of 
(13.581) for = /3fe = 1 by Z\ and for bk = j5k = —f hy Z 2 . By replacing Vk^ Sk, bk and 
/3fc by suitable Vk, Sk, bk and f3k respectively, we may assume that 

(3.63) zi= Xl + (heskVk + f{sk)vkX^) 
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with |sfc| < where 0 < /i < 1. We note that if we change the signs of bk and Pk on 
the RHS of (I3.63P then this RHS is equal to Z 2 . We explain why p.63l) holds. Let us 
denote for 0 < a < b the following piece of a graph by 


(3.64) r(a, b) = {(es, e^/(s)) : a < s < b} . 

The minimum of dist(0, •) on T = r(e“b e““^) is attained in (es, e^/(s)) where 


(3.65) 


s = 



here and in the following we denote by the origin 0 the point X^. Furthermore, the 
distance of the left endpoint of F to the origin is larger than the distance of the right 
endpoint of F to the origin, to see this we omit constants and check the orders of e in 
the squares of these distances. We get 

(3.66) g2+2ai 


for the left endpoint and 

(3.67) ^2-2a2 ^ 


for the right endpoint. Therefore we must have 


(3.68) 


4- 


ai 


4j 

1-7 


<2 - 2a2 


which follows from (12.91) . Hence all points p £ r(s, e““^) can be ’covered’ by rotating 
F(e“i,s) around the origin, i.e. there is a rotation Op of the plane around the origin 
(depending on p) such that 


(3.69) pe CipF(e“b5). 

But this rotation can be realized by choosing the quantities with a tilde suitable. 

Now, we look at the RHS of (|3.61l) in the situation that Sk,Vk,bk, Pk carry a tilde. 
The only possibly ’bad’ summands therein are the summand which contains bk and 
the one with Pk (note, that Sk is small). Let us denote the first by A and the second 
by H. If H H > 0 we set and if H H < 0 we set X^j^^ = Z 2 - 

(ii) Let us consider the case where X(j){xo,to) ^ 0 and set do = \V(j){xo,to)\- We 
let p, —p, q and —q be the four possible points for X^^^ depending on the choice of 
bk,Pk where we let p be the point with bk = Pk = ^- We choose bk,Pk so that the 
summands in which they appear are non-negative. Let cq > s be a sufficiently small 
fixed constant (not depending on e) then (I3.62|) follows from (13.611) as can be seen 
from the following cases. 

Case (1): We assume Sk < cq. Let us denote the summands of the resulting 
RHS of (13.611) from left to right by H, B, C, and D. If | (V(p(2f^), | > ^ then 

A > is the dominating summand, otherwise | {Vip{Xl),Vk) | > ^ so 

that C > OoP^ dominates B and D and hence the claim follows. 

Case (2): We assume Sk > cq. 

Case (2a): Assume that there is a sufficiently large constant ci > 0 so that 

(XpiXD^vi) > Ciesfe. 


(3.70) 
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Then the first summand on the RHS of (13.611) is dominating and gives the desired 
estimate. 

Case (2b): Let 


(3.71) 


< ciesfc- 


Case (2ba): Let curv((/))|x' < 0. 
We have 


(3.72) 


vi = 


|v</>(w^)| 


+ cesfc 


where we possibly have to change the signs of both bk and /3fc. We denote the sum¬ 
mands of the RHS of (13.611) from left to right by B, C and D. We have 

(3.73) |H| < cehl 


and 


(3.74) 


e-\B + C + D) 

>dtc)ixi) + fis,)\y4>ix^)\ 





\v<j^{xi)y \vc^{xi)\ / 


-cesl 


=dtHXl) + fis,)\y4>iX^)\- |V^(X^)|| curv(<^)|;f. - cesl 
=dt<l){Xk) - |V0(X^)|v?«,(curv(^i)|x^) - cesl 
>dt<l>{Xl) - |V(/)(X^)|(/?(curv((/))|x») - cesl 


Now, the claim follows from assumption (12.91) . 

Case (2bb): Let curv(()))|jfe > 0. W.l.o.g. we may assume that 

(3-75) curv(()))|(,j,„_iQ) > 0, 


otherwise we choose a smaller ^-neighborhood. Then dt4>{Xy) > We define 

(3.76) Vk = OpVk, vy = OpV^ 

where Op as in (13.691) (and p the point corresponding to = /3fc = 1) then 

(3.77) p = Xl + (eSfeUfc -b e^f{sk)vk,e'^) 


with suitable small Sfc. Then 

(3.78) |(V<^(X^),u^)| >co>0. 

We denote the terms on the RHS of (13.611) (now in the situation of the quantities with 
a tilde) from left to right by A, B, C and D. There holds > |I1|. By considering 
—p instead of p we may assume that H -|- C > 0. Since B > e^^ the claim follows. 
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We used that the rotation angle a of Op is almost hence we supplement an 
estimate for a. We write a = a 2 — cti with 


(3.79) 
and 

(3.80) 

where 

(3.81) 

or, equivalently, 

(3.82) 
so that 

(3.83) 

Altogether this gives 

(3.84) 


tanai = 


eV(sfc) 




= ces.: 


k ’ 


Sk > Co 


tanQ;2 = 


V(gfc) 

esfc 


Sk < s 


e^/(sfe) > cesfc, 


Sk < ce 27 s 


tana2 > e s^'’ 


1+^ 

«i = 0(esj:-") A 


a2- 2 


1-7 _i±Tl 
= 0 ( e—Sfc "" 


4. Numerical scheme. Let e > 0 be a small step size and h > 0 the spatial 
step size. We consider the points flh = hi? forming a rectangular grid in the plane. 
Our value function is given at the final time T > 0, 


(4.1) 


‘^i;T)=uo 


where we assume that uq = 0 in \ 0. The calculation of the value function is 
backward in time. Let x S flh be a grid point. In order to calculate u'^{x, t — e^) from 
we propose the following stategy. 

We discretize the control set A by 


(4.2) 


jAs ^ I^ai + ^As : 0 < r < m}. As = 


g—a2 _ gOl 

ro 


and the set of unit vectors {r;GM^:||r!|| = l}by the discrete subset 


(4.3) 


QiTtI 

Sig = {(cos a, sin a)) : a = ——,0 < Z < ^o} 

h 


where rg, Iq are sufficiently large natural numbers. We then set 


(4.4) 


u^ix^t — e^) = min max u^{z,t) 


z = X + besv + /3e^/(s)u'* 


where 

(4.5) 
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Table 1: Error for e = 0.08 and 7 = 0.7 


scale = 0.1 

scale = 0.3 

scale = 0.9 

loo-error 0.1085 

0.08432 

0.08431 

Zi-error 0.1571 

0.1006 

0.0988 


and u'^ on the RHS of (14.41) is replaced by the value of a bilinear interpolating function 
if z is not a grid point. 

Let us look at the relative scaling of e and h in the special case of 7 = 0.5. From 
(13.651) and (13.661) we know that the distance from z to a; is at least ces and at most 
ce^(^-“i), note, that ^ < Oi < 1 and 0 < 02 < min(2ai — 1, i). In order to include 
second order information of the value function in the optimization problem (14.41) it 
is natural to assume that the range of the control variable z is not covered by a 

o 4 

square of size h , so that we assume h < e 3 . Then we get for the time step size the 
estimate and the standard interpolation error in the supremum norm is of 

order = o(e^). 

5. Numerical examples. We consider the case of a shrinking circle with initial 
radius Rq = 1 and different values of | < 7 < 1. As initial function we use 

(5.1) uo(x) = max(|a;|'’''''^ — 
with solution 

(5.2) u{x, t) = max(|a;|'’''''^ — + (7 + 1)^, 0)^ 

for 0 < t < Tmax = To define the value function completely we have to 

specify the parameters e,ai and 02 • For the numerical approximation of we have 
to specify the discretization parameters h,lo,rQ. 

5.1. Search of good parameters 01,02 for the example case e = 0.08 and 

7 = 0.7. In this subsection we consider the case e = 0.08 and 7 = 0.7 (then T^ax ~ 
0.59) and test different values of 01,02 while the numerical parameters h = 0.01, 
Iq = 160 are fixed and tq is adapted so that A is discretized equidistantly with step 
size 0 . 01 . 

In order to determine good values for oi and 02 we denote the lower bound for 
oi in (12.81) by rui, the upper bound by m 2 and the upper bound for 02 in (j2.9ll by 
m 3 . We introduce a scaling factor scale and set 

(5.3) oi = mi + scale •(m 2 — mi) A 02 = scale-ma. 

The approximation error, i.e. the difference u — u'^, is measured with respect to 
the Zoo(hZ^)-norm and with respect to the Zi(/iZ^)-norm (more precisely, we scale the 
latter norm by h?) at each time step. Then the supremum over all time steps until 
time T = 0.12 of both norms is taken and denoted by ^oo-error and /i-error. 

Table [T] conhrms the expected behavior that a larger control interval A leads to 
better convergence properties. 
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Table 2: Influence of h on the error for e = 0.08, 7 = 0.7, scale = 0.9, Iq = 160 and ro = 160 




/ 00 -error 

/i-error 

h 

= 0.16 

0.2639 

0.1876 

h 

= 0.08 

0.1180 

0.0999 

h 

= 0.04 

0.0947 

0.1035 

h 

= 0.02 

0.0868 

0.1045 


Table 3: Influence of ro on the error for e = 0.08, 7 = 0.7, scale = 0.9, h — 0.01 and 
lo = 160 


/ 00 -error /i-error 


ro 

= 10 

0.0902 

0.1075 

ro 

= 20 

0.0860 

0.1014 

ro 

= 40 

0.0848 

0.0998 

ro 

= 80 

0.0844 

0.0991 

ro 

= 160 

0.0843 

0.0988 


5.2. Influence of the numerical parameters h, Iq and rg. We consider for 
the case e = 0.08, 7 = 0.7 and scale = 0.9 three scenarios in which we analyze the 
influence of the numerical parameters h, rg and /g on the accuracy of our algorithm. 
Therefore we keep hold two of these three numerical parameters in each scenario 
while the remaining one runs through some test values. The corresponding errors are 
presented in Tables [H |3] and [H 

5.3. Convergence of the value functions. In Tables [5] and [H] we present the 
approximation errors in dependence of e for 7 = 0.8 and 7 = 0.9 (we set scale = 0.9). 

5.4. Zero level sets of the value function. Figure [T] shows the zero level sets 
of u'^ for different times. The parameters are chosen as follows: 7 = 0.9, e = 0.04, 
h = 0.02, scale = 0.9, rg = 100 and /g = 160. Figure [5] shows the same but now the 
initial function is replaced by 

(5.4) {tg(a;) = max ^(xi -I- 1 . 7 a; 2 )^ — 1,0^ , x=[xi,X 2 ) 

so that the initial curve is an ellipse. One can see that the evolving curves, i.e. the 
zero level sets of u'^, become circular and shrink to a point which is the known behavior 
for the level sets of the limit function u. 

5.5. Discussion of the numerical results. We recall that there are four nu¬ 

merical parameters, i.e. h, e, rg and /g, in our approach. Incorporating the two 
space dimensions we are confronted with a de facto flve dimensional problem which 
we discretize. Still, our best error is 0.0139 and in Tabled attained for the moderate 
parameter values h = 0.01, At = = 0.0004, rg = 100 and /g = 360. We observed 

that the convergence properties improve when 7 is closer to 1 . 

To classify our errors we compare with m where for the Kohn-Serfaty scheme 
of the curve shortening flow a / 00 -error of 0.0078 is obtained where At = 0.02 and 
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Table 4: Influence of Ig on the error for e = 0.08, 7 = 0.7, scale = 0.9, h = 0.01 and 
ro = 160 




/oo-error 

/i-error 

^0 

= 10 

0.4254 

0.4333 

^0 

= 20 

0.1885 

0.2463 

^0 

= 40 

0.1291 

0.1327 

lo 

= 80 

0.0922 

0.1114 

lo 

= 160 

0.0843 

0.0988 


Table 5: Convergence of the value functions m'' for 7 = 0.8 



/oo-error 

/i-error 

h 

ro 

^0 

e = 0.09 

0.0339 

0.0261 

0.01 

100 

360 

e = 0.08 

0.0325 

0.0253 

0.01 

100 

360 

e = 0.05 

0.0250 

0.0170 

0.01 

100 

360 

e = 0.04 

0.0205 

0.0112 

0.01 

100 

360 

e = 0.02 

0.0139 

0.0130 

0.01 

100 

360 


Ax = 0.0098 , cf. [11] Table 2]. For the same flow the semi-Lagrangian scheme leads 
to a loo“Crror of 0.00832, cf. m Table 2]. Note, that therein the initial function is a 
higher order power of the difference \x\^ — Rq than in our paper and this effects the 
error. Furthermore, we refer to m Table 1] where a semi-Lagrangian approximation 
for our equation (11.21) in the case 7 = | is considered. For Ax = 0.01 and At = 0.001 
the authors obtain there a / 00 -error of 0.0542 and of 0.0306 for a modified scheme. 
Note, that the time interval for the calculation in each of the above cases is—as in our 
case—approximately 20% of the time interval [0,T*) where T* is the time at which 
the evolving curve ’becomes a point’. Hence our Tables [5] and [6| contain values of 
adequate size. 
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